098 


DDC 

■j^rarnnoji 

J  JUL2  4I983 

JtifhOED  \J  [ 


TISIA  B 


LMSC-A324768 
10  June  1963 


PROGRAM  461  TELEMETRY 
ALIASING  ERROR  ANALYSIS 


APPROVED: 

y  .  / 

E.  SCHAEFE R",  MANAGE R 
PROGRAM  461  C&C  SYSTEM 


APPROVED: 


K 


/  y 

/  C. SOLVASON,  MANAGER 
PROGRAM  461 


LOCKHEED  MISSILES  8c  SPACE  COMPANY 


LMSC-A.324768 


FOREWORD 


This  report  is  published  by  Lockheed  Missiles  and  Space 
Company,  Space  Systems  Division,  under  Contract 
AF  04(647) -787 .  It  presents  the  results  of  a  computer 
simulation  of  Program  461  payload  telemetry  system  in  a 
form  which  enables  the  designer  to  determine  data  sampling 
rate  and  number  of  channels  for  a  given  aliasing  error 
allocation. 
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INTRODUCTION  AND  SUMMARY 


Aliasing  error  occurs  when  interpolation  is  applied  to  data  samples  whose  period 
between  samples  is  too  great  to  recover  the  higher  frequency  components  in  the  data 
spectrum.  Since  the  data  spectrum  generally  contains  high  frequency  components 
which  are  not  of  interest,  the  signal  is  low-pass  filtered  to  attenuate  the  undesirable 
portion  of  the  spectrum.  It  is  this  portion,  although  attenuated,  that  "folds"  about  the 
Nyquist  frequency  thereby  causing  aliasing  error. 

Aliasing  error  may  be  reduced  by  faster  sampling  and/or  better  filtering.  Since  it  is 
mandatory  that  a  high  degree  of  accuracy  be  attained  with  a  minimum  sampling  rate, 
an  optimization  analysis  is  clearly  necessary.  To  accomplish  this,  the  telemetry 
system  is  described  mathematically  by  the  transfer  functions  of  its  significant  blocks . 
Previously,  a  cursory  analysis  was  performed  by  using  graphical  integration  (Ref  1) . 
This  effort  led  to  a  more  mathematically  refined  approach  with  the  calculations  being 
performed  by  an  IBM  7094  computer. 

Mathematical  functions  are  derived  for  the  Program  461  payload  telemetry  system, 
and  utilized  to  compute  rms  aliasing  error  versus  data  sampling  rate.  Three  dif¬ 
ferent  interpolation  filter  types  are  compared;  a  physically  realizable  Butterworth, 
an  ideal  rectangular  shape  with  infinite  slope,  and  an  optimum  filter  whose  charac¬ 
teristics  are  a  function  of  the  data  spectrum.  Both  the  high  resolution  and  low  resolu¬ 
tion  channels  are  analyzed.  The  computation  was  performed  on  an  IBM  7094  computer, 
with  a  program  written  so  that  system  parameters  could  be  easily  varied. 

For  the  two  blur-circle  time  functions  assumed,  Gaussian  and  cosine-squared,  curves 
of  rms  aliasing  error  versus  sampling  frequency  are  developed.  Finally,  for  a  given 
aliasing  error,  interpolation  filter,  and  input  spectrum,  the  allowable  number  of  low 
and  high  resolution  channels  are  shown. 
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SECTION  1 

DESCRIPTION  OF  TELEMETRY  SYSTEM 

1.1  BLOCK  DIAGRAM 

Figure  1-1  shows  the  significant  blocks  of  the  system ,  the  voltage  transfer  functions , 
and  the  sampling  functions . 

1.2  OPTICS 

When  a  point  source  is  viewed  by  the  optical  system,  the  energy  distribution  of  the 
image  forms  a  blur-circle  whose  mathematical  expression  is  a  function  of  the  optical 
parameters.  As  the  blur-circle  moves  across  the  detector  aperture,  the  energy 
versus  time  curve  assumes  a  shape  somewhere  between  a  Gaussian  and  a  cosine- 
squared  pulse.  Figure  1-2  shows  the  frequency  spectrum  for  both  cases. 

1.2.1  Case  I  -  Gaussian  Pulse 


h 


o 


-  co<  t  <  +  CO 


It  was  assumed  (Ref  1)  that  since  95  percent  of  the  image  is  contained  within  the 

1  /I 

detector  aperture  and  the  scan  time  across  the  aperture  is  -^q  seconds,  t^  =  1/4  (j2Q 
This  is  due  to  the  fact  that  95  percent  of  the  energy  in  a  Gaussian  pulse  is  contained 
within  -2  a  <t  <  +2  <r.  The  frequency  function  or  spectrum  of  the  Gaussian  pulse  is 
also  Gaussian  and  equal  to: 
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Figure  1-1  System  Block  Diagram 
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Figure  1-2  Frequency  Spectrum  of  Assumed  Blur-Circle 
Time  Functions 
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(f)  = 


1.2.2  Case  II  -  Cosine-Squared-Pulse 


>0  (0  = 


2A  2  t 

—  C°S  7T  —  , 


t  t 

<  t  <  +  — 
2  ~  2 


=  0 ,  elsewhere 

Again  for  this  case,  95  percent  of  the  image  is  assumed  to  be  contained  within  the 
detector  aperture;  consequently  tQ  =  (Appendix  A) .  The  normalized  frequency 
function  is: 


490  A 


.  TTf 
Sm  490 


r  -  xy  1 

L  \490 / 


1.3  DETECTOR  APERTURE 


Since  the  detector  shape  is  rectangular,  a  sweep  across  the  aperture  would  result 
in  a  rectangular  time  function: 


M-f'  -f 


=  0,  elsewhere 
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Since  the  sweep  across  the  aperture  is  seconds  for  the  high  resolution  channels 

1  ^11 
and  —  seconds  for  the  low  resolution  channels,  t  =  and  ^7:  respectively.  The 
obu  0  i  obU 

frequency  function  is: 


ha  (f)  "  A 

t 

O 

1.4  DETECTOR  AND  AMPLIFIER 

The  frequency  function,  in  the  case  of  the  high  resolution  channels,  for  the  detector- 
amplifier  combination  was  derived  from  the  idealistic  shape  as  shown  in  Figure  1-3. 
Two  situations  were  postulated;  one  in  which  the  high  frequency  roll-off  of  the  existing 
amplifier  is  18  db/octave,  and  one  in  which  the  roll-off  is  increased  to  24  db/octave. 
Since  the  amplifier  is  composed  of  R-C  coupling  networks ,  the  equation  for  the  gain 
across  the  band  is  easily  formed.  For  N  =  3  and  N  =  4,  the  high  frequency  -3  db 
point  is  214  cps  and  183  cps  respectively.  However,  the  curve  is  still  asymptotic  to 
the  slope  as  illustrated  in  Figure  1-2. 
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Where  N  =  3  for  18  db/octave  roll-off,  and  N  =  4  for  24  db/octave  roll-off. 

In  the  case  of  the  low  resolution  channels,  the  high  frequency  pole  is  reduced  from 
420  to  210  cps. 


1.5  MULTIPLEXER 


Ideally,  the  multiplexer  samples  the  data  repetitively  with  zero  dwell  time  thereby 
producing  a  series  of  delta  functions  of  magnitude  equal  to  the  magnitude  of  the  data 
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at  the  sampling  time.  Consequently,  if  only  the  peak  value  of  the  data  is  desired,  the 
error  would  be  zero  if  the  effect  of  finite  dwell  time  is  not  considered.  However,  since 
the  spectrum  of  the  data  is  to  be  reproduced,  interpolation  between  data  samples  must 
be  applied  and  error  is  inherent.  The  effect  of  finite  dwell  time  is  to  introduce  fre¬ 
quency  distortion  in  the  recovered  data.  The  method  of  determining  the  magnitude  of 
the  error  is  explained  in  Reference  1.  The  error  is  0.004  percent,  with  typical  system 
parameters  (Appendix  B) .  The  sampling  function  ideally  is: 

K(<)-  2  Ss(f-nfJ 

n  =0 


where 

Sg  (t'j  =  Sampled  data  power  system 

Sg  ^  f  J  =  data  power  spectrum 

f  =  sampling  frequency 
s 

1.6  DEMULTIPLEXER 

The  demultiplexer  merely  decommutates  the  composite  signal  and  has  no  effect  upon 
the  aliasing  error . 

1.7  INTERPOLATION  FILTER 

The  three  interpolation  filter  characteristics,  which  were  considered,  are  listed 
below. 

1.7.1  Ideal  Filter 

The  ideal,  rectangular  shaped  filter,  with  infinite  slope  at  half  the  sampling  frequency, 
is  easy  to  simulate  with  the  computer.  Its  function  is: 
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H 


o 


f 

1,  0  <  f  <  -f 


=  0 ,  elsewhere 


1.7.2  Butterworth  Filter 

A  10-pole  60  db/octave  Butterworth  filter,  which  is  physically  realizable,  is  given  by 
the  function: 


1.7.3  Optimum  Filter 


The  optimum  filter,  as  derived  by  Spilker  (Ref  2),  is  of  the  form: 


Realization  of  this  filter  depends  upon  a-priori  knowledge  of  the  data  spectrum;  con¬ 
sequently,  we  may  only  approach  optimum  in  actual  practice. 
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SECTION  2 

METHOD  OF  ANALYSIS 


The  method  of  analysis  following  that  of  J.  J.  Downing's  work  (Ref  3),  along  with  the 
system  transfer  functions  described  previously,  gives  an  exact  calculation  of  the 
aliasing  error . 

2.1  MEAN-SQUARED  ERROR 


The  sampled  data  spectrum  may  be  written  as: 

(f)  -  Ss  (f)  +  Sn  (0 

where 

CO 

\(l)  -  E  Ss(f-"fs) 
n=  1 


For  a  general  interpolation  filter  vnth  voltage  transfer  function  Hc  ^f  J  ,  the  mean- 
squared  error  is: 


For  the  optimum  interpolation  filter,  the  mean-squared  error  is: 
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where  the  optimum  filter  transfer  function  is  given  by: 


H 

c 


2.2  NORMALIZED  ERROR 


The  error  was  normalized  with  respect  to  the  data  power  spectrum 


and  expressed  as  percent  rms  aliasing  error  in  the  subsequent  curves. 
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SECTION  3 

RESULTS  OF  COMPUTER  ANALYSIS 


To  facilitate  IBM  7094  computer  analysis  the  sampled  data  spectrum  was  summed  to 

2.5  f  and  the  mean-squared  error  was  integrated  from  0  to  2.5  f  and  multiplied  by 

2  as  allowed  by  symmetry.  Since  the  roll-off  of  the  assumed  interpolation  filter  is 

quite  sharp  (infinite  in  the  case  of  the  ideal  filter) ,  negligible  improvement  is  gained 

by  summing  and  integrating  further  than  2.5  f  .  The  operations  were  performed  in 

s 

50  cycle  increments  in  order  to  give  realistic  results  while  keeping  machine  time 
reasonable . 

3.1  ALIASING  ERROR  VS  SAMPLING  FREQUENCY 

3.1.1  High  Resolution  Channels 

Figures  3-1  and  3-2  show  the  rms  aliasing  error  for  three  interpolation  filters  when 
the  input  time  function  is  Gaussian.  Figures  3-3  and  3-4  show  the  rms  aliasing  error 
for  the  three  interpolation  filters  when  the  input  time  function  is  cosine-squared. 

The  rms  aliasing  error  is  slightly  less  for  the  cosine-squared  function,  due  to  less 
high  frequency  energy  in  the  data  spectrum,  than  for  the  Gaussian  function.  This 
comparison  may  be  made  by  referring  to  Figure  1-2,  and  noting  the  difference  between 
the  two  curves . 

3.1.2  Low  Resolution  Channels 

Figures  3-5  and  3-6  show  the  results  of  the  analysis  when  the  aperture  and  amplifier 
functions  are  modified  for  the  case  of  low  resolution  channels  (Gaussian  input  pulse) . 
Figures  3-7  and  3-8  show  the  rms  aliasing  error  for  the  low  resolution  channels  when 
the  input  time  function  is  cosine-squared. 
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3.2  EFFECT  OF  CHANGING  AMPLIFIER  ROLL-OFF 

If  the  amplifier-detector  high  frequency  roll-off  is  increased  to  24  db/ octave,  a 
significant  reduction  of  error  is  achieved.  Figure  3-9  shows  the  results  of  the 
computer  analysis  for  this  assumption  with  a  Gaussian  input  function  to  the  high 
resolution  channels,  and  an  ideal  interpolation  filter.  For  the  same  sampling  fre¬ 
quency,  the  aliasing  error  is  reduced  approximately  50  percent.  This  should  have 
been  immediately  obvious  since  the  area  under  this  portion  of  the  curve  is  reduced 
by  one-half  when  the  slope  is  doubled. 

3.3  EFFECT  OF  CHANGING  THE  INTERPOLATION  FILTER  CUT-OFF 

For  a  constant  sampling  frequency,  the  cut-off  frequency  of  the  ideal  and  Butterworth 
filters  was  varied  about  half  the  sampling  frequency.  The  computer  results  indicated 
that  minimum  error  occurred  for  a  cut-off  at  fgy{, . 
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Figure  3-9  Effect  of  Increasing  Amplifier-Detector  Roll-Off  (Aliasing 
Error  vs  Sampling  Frequency  -  High  Resolution  Channels: 
Gaussian  Input  Function) 
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SECTION  4 

NUMBER  OF  HIGH  VS  NUMBER  OF 
LOW  RESOLUTION  CHANNELS 


The  data  spectra  used  in  this  study  are  typical  of  the  high  and  low  resolution  sensors, 
and  since  high  accuracy  is  desired,  the  aliasing  error  must  be  minimized.  Since  the 
Gaussian  input  function  is  a  "worst-case"  condition,  it  will  be  used  in  calculating  the 
number  of  channels. 

For  a  seven-bit  PCM  system,  the  maximum  commutation  rate  is  114  Kc  as  dictated 
by  the  800  kilo-bit  maximum  rate  set  by  IRIG.  Therefore,  the  number  of  channels 
that  can  be  multiplexed  is: 

xt  _  114,000 


where  f  is  the  sampling  rate  per  channel . 
s 

Since  f  is  dictated  by  the  tolerable  aliasing  error,  the  number  of  channels  could 
s 

easily  be  determined  were  it  not  for  the  difference  in  sampling  rates  required  by  the 
high  and  low  resolution  channels  for  the  same  aliasing  error. 

For  example:  if  0.3  percent  rms  aliasing  error  is  tolerable  it  is  indicated  (Figure  3-2) 
that  1700  sps  should  be  used  for  the  high  resolution  channels,  and  that  1190  sps  (Fig¬ 
ure  3-5)  should  be  used  for  the  low  resolution  channels.  Therefore,  in  order  to 
operate  in  an  efficient  manner,  the  high  and  low  resolution  channels  should  be  sampled 
at  a  rate  of  or  y.  Because  this  would  require  an  overly  complex  multiplexer 
design,  the  1:1  sampling  ratio  (although  inefficient)  is  desirable.  With  this  arrange¬ 
ment,  the  number  of  channels  is  dictated  by  the  tolerable  aliasing  error  of  the  high 
resolution  channels  For  a  0.3  percent  error  the  total  number  of  channels  is  67 . 
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Super-commutation  may  be  employed  providing  the  allowable  aliasing  error  for  the 
low  resolution  channels  is  relaxed  to  that  value  corresponding  to  half  the  sampling 
rate  of  the  high  resolution  channels.  If  0.3  percent  aliasing  error  is  tolerable,  the 
sampling  rate  (Figure  3-2)  is  1700;  therefore,  at  a  sampling  rate  of  850  for  the  low 
resolution  channels  the  aliasing  error  for  the  channels  (Figure  3-5)  is  1.7  percent. 

If  super -commutation,  twice-per-frame ,  is  utilized  for  the  high  resolution  channels, 
the  total  number  of  channels  is  given  by 


N  =  2H  +  L 

Where  H  is  the  number  of  high  resolution  channels,  N  is  the  number  of  low  resolution 
channels,  and  the  total  number  of  channels  is  determined  by  the  sampling  rate  of  the 
low  frequency  channels.  In  the  case  cited: 


N 


114,000 

850 


134 


Figure  4-1  indicates  the  number  of  high  vs  the  number  of  low  resolution  channels  for 
the  foregoing  assumption. 
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SECTION  5 

CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  CONCLUSIONS 

Results  of  the  analysis  indicate  that  the  assumption  of  a  Gaussian  function  for  the  blur- 
circle  gives  a  "worst-case"  condition.  Therefore,  pertinent  conclusions  should  be 
drawn  from  Figures  3-2,  3-5,  and  3-6  using  the  curve  for  the  ideal  interpolation  filter. 
It  is  also  indicated  that  super-commutation,  twice-per-frame,  can  be  employed  if  the 
accuracy  requirements  are  relaxed  for  the  low  resolution  channels  thereby  allowing  a 
greater  total  number  of  channels.  If  once -per -frame  commutation  is  used,  the  total 
number  of  channels  is  67  for  a  0.3  percent  maximum  aliasing  error  at  114  Kc  com¬ 
mutation  rate .  This  condition  is  wasteful  of  bandwidth  if  utilized  to  telemeter  the  low 
resolution  channels.  When  twice-per-frame  super-commutation  is  employed,  the 
total  number  of  channels  is  increased  to  134,  but  the  aliasing  error  of  the  low  resolu¬ 
tion  channels  is  increased  to  1.7  percent. 

By  modifying  the  transfer  function  of  the  detector-amplifier  to  increase  the  roll-off  to 
24  db/octave,  the  aliasing  error  is  reduced  approximately  50  percent. 

5.2  RECOMMENDATIONS 

Since  the  study  is  based  on  the  assumption  that  the  actual  blur  circle  function  is  some¬ 
where  between  Gaussian  and  cosine-squared,  it  is  recommended  that  a  laboratory 
experiment  be  devised  to  determine  the  actual  function.  The  amplifier-detector  trans¬ 
fer  function  could  also  be  made  more  precise  if  the  unit  were  constructed  and  frequency 
response  measurements  performed.  With  these  improvements  in  the  transfer  functions 
the  aliasing  error  analysis  should  be  re-run  with  other  refinements  such  as  inclusion 
of  detector  and  thermal  noise . 
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Under  consideration  at  the  present  time  is  the  effect  of  detector  noise  on  the  results 
of  the  analysis.  Work  should  continue  in  order  to  determine  the  mathematical  function 
which  best  describes  the  noise  spectral  density  and  how  to  treat  the  combination  of 
signal  and  noise  statistically  as  it  is  passed  through  the  system. 

Consideration  should  also  be  given  to  increasing  the  roll-off  of  the  amplifiers  from 
18  db/octave  to  24  or  30  db/octave.  It  is  recognized  that  the  penalty  for  doing  this  is 
increased  complexity  with  all  its  ramifications,  but  the  outcome  is  great,  as  witnessed 
by  the  curves  in  Figure  3-5. 

Since  the  total  number  of  telemetry  channels  is  sharply  reduced  if  the  accuracy  re¬ 
quirements  are  identical  for  the  high  and  low  resolution  channels,  the  feasibility  of 
relaxing  the  low  resolution  accuracy  requirement  should  be  investigated. 
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Given: 


APPENDIX  A 

DERIVATION  OF  THE  PARAMETERS 
FOR  COSINE-SQUARED  PULSE 


=  0,  elsewhere 


95  percent  of  energy  (area)  lies  between  -b  and  +b,  and  the  time  between 
-b  and  +b  =  -=7-7:  second 


Find: 


t  in  terms  of  b  ^ 

o 


Solution: 


Total  Area  =  A 


A 

2 


t 

_o 

2 


A 


t 

0 

2 
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Area  between  -b  and  +b  = 


A 

2 


,  ,  o  .  2tt 

t+2W  Sln  t-* 
o 


+b 


-b 


0.95  Ay 


=  A 


,  ,  o  .  2tt 

t+  2?f  Sin— 4 

o  J 


0.95  Ay  =  A 


,  *o  .  2rr  , 

b  +  2?  Sm  t"  b 
o 


This  is  a  transcendental  function;  for  which  the  approximate  solution  is: 


b  =  0.34  t 

o 


Since  the  time  between  -b  and  +b  is  yy  second; 


t 

o 


1 

490 


second 
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APPENDIX  B 

EFFECT  OF  FINITE  DWELL  TIME 


It  is  shown  (Ref  1)  that  the  error  caused  by  frequency  distortion,  introduced  by  finite 
sampling,  is  given  by: 


sin  tt  f  t 
TT  f  T 


Where: 

f  is  the  highest  frequency  of  interest  in  the  data,  which  in  this  case  is  assumed 
to  be  the  Nyquist  frequency  fgy^,  typically  850  cps;  and  t  is  a  sample  duration. 
For  a  114  Kc  commutation  rate,  100  percent  duty  cycle  multiplexer: 

T  =  -  1  ■ 

114,000 

Therefore,  the  error  at  1840  cps  is: 

1  -  0.99996  or  0.004  percent 
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